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ABSTRACT 

A variety of on-orbit imaging and spectroscopic observations are used to char- 
acterize the Charge Transfer Efficiency (CTE) of the Charge- Coupled Device 
(CCD) of the Space Telescope Imaging Spectrograph (STIS) aboard the Hubble 
Space Telescope. A set of formulae is presented to correct observations of point 
sources for CTE-related loss of signal. For data taken in imaging mode, the CTE 
loss is parametrized in terms of the location of the source on the CCD, the source 
signal level within the measurement aperture, the background level, and the time 
of observation. For spectroscopic data, it is found that one additional parameter 
is needed to provide an adequate calibration of the CTE loss, namely the signal 
in the point spread function located between the signal extraction box and the 
read-out amplifier. The effect of the latter parameter is significant for spectra 
taken using the G750L or G750M gratings of STIS. The algorithms presented 
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here correct flux calibration inaccuracies due to CTE losses as large as 30% to 
within ~ 1.5% RMS throughout the wavelength range covered by the STIS CCD 
modes. This uncertainty is similar to the Poisson noise associated with a source 
detected at a signal level of about 2500 electrons per resolution element. Using 
bi-directional CCD readouts, centroid shifts incurred due to CTE loss are also 
derived. A tight correlation is found between the CTE loss and the centroid 
shift (both for imaging and spectroscopic modes), thus enabling one to correct 
for both effects of imperfect charge transfer to STIS CCD observations. 

Subject headings: instrumentation: detectors — methods: data analysis 



Introduction 



Astronomical observation was revolutionized about 25 years ago by charge-coupled de- 
vice (CCD) technology, due to a combination of generally linear response over a very large 
dynamic range and high quantum efficiency However, a shortcoming of CCDs is the fact 
that the transfer of charge from one pixel to the next is not perfect. Charge Transfer Ef- 
ficiency (CTE) is the term commonly used to describe this effect, and it is quantified by 
the fraction of charge successfully moved (clocked) between adjacent pixels. In practice it is 
often more useful to use the term Charge Transfer Inefficiency (CTI = 1— CTE). We will do 
so in this paper when considering quantitative measures. The main observational effect of 
CTE loss is that a star whose induced charge has to traverse many pixels before being read 
out appears to be fainter than the same star observed near the read-out amplifier. 

Laboratory tests have shown that the CTE loss of CCDs in creases significantly when 
being subjected to radiation damage (e.g., IJanesick et al.lll99lt ). This is particularly rel- 
evant for space-borne CCDs such as those aboard Hubble Space Telescope (HST), where 
the flux of high-energy particles, particularly in the South Atlantic Anomaly, is signifi- 
cantly higher than on the ground. The effect is significant for all CCD detectors used 
on HST instruments (e.g., WFPC2: Whitmore, Heyer, & Casertano 1999; Dolphin 2000; 
STIS: Goudfrooij & Kimble 2003; ACS: Riess & Mack 2004), and charge losses increase 
with elapsed on-orbit time. The purpose of the current paper is to characterize the CTE 
loss of the CCD of the Space Telescope Imaging Spectrograph (STIS) for observations of 
point sources in terms of its dependencies on the X and Y positions, target intensity, back- 
ground counts, measurement aperture size, observing mode, and elapsed on-orbit time. Sev- 
eral aspec ts of on-orbit characterizations o f the CTE loss of the STIS CCD have been re - 
ported bv iGilliland. Goudfrooii. fc Kirnbl] Jl999h. iKimble. Goudfrooii. fc Gillilandl (boooh 
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Goudfrooij fc Kimbld (120031 ) , and iBohlin fc Goudfrooijl (120031 ) . The current paper provides 
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a comprehensive description of the final STIS CTE calibrations, using data taken through 
mid-2004 (i.e., just before STIS operations ended) which allows a more accurate temporal 
dependence of the CTE loss than the reports mentioned above. The current paper also 
includes analysis of data taken specifically to constrain the CTE loss at lower signal levels 
than before. Furthermore, we provide (for the first time) (i) separate characterizations of 
the CTE loss for spectroscopic and imaging modes, (ii) a comparison between data taken 
in gain = 1 and gain = 4 e~ DN _1 settings, and (Hi) the relation of the CTE loss with its 
associated centroid shifts. 

The STIS CCD is a 1024 x 1024 pixel, backside-illuminated device with 21 /mi x 21 /mi 
pixels. The CCD was built by Scientific Imaging Technology (SITe) with a coating process 
that allows coverage of the 200-1000 nm wavelength range for STIS in a wide variety of 
imaging and spectroscopic modes. Key features of the STIS CCD architecture are shown 
schematically in Figure [TJ Two serial registers are available. A read-out amplifier is located 
at all four corners, each with an independent analog signal processing chain. The full image 
can be read out through any one of the four amplifiers, or through two- and four-amplifier 
combinations. By default, science exposures employ full- frame readout through amplifier 
'D', which features the lowest read-out noise (4.0 e~ RMS at launch). 



Further technical details regarding the STIS CCD in particular are provided in lKimble et al 



( 1994). while backgrou nd information on the design of STIS in general can be found in 



Woodgate et all (Il998h . 



This paper is organized as follows. Section [2] describes methods used to measure the 
CTE loss: Two methods designed to quantify different aspects of CTE loss for observations 
of point sources in sparse fields for spectroscopic and imaging modes, and one method to 
improve the characterization of the CTE loss for point source spectroscopy, using on-orbit 
observations of spectrophotometric standard stars. Section [3] presents the results of the CTE 
analyses, and Section [4] summarizes the results. 



2. CTE Measurement Methods 
2.1. Spectroscopic Modes 

2.1.1. Internal Sparse Field Test 

This novel test method, designated the "internal sparse field" test, was developed by the 
STIS Instrument Definition Team. The method quantifies two key aspects of CTE effects 
on spectroscopic measurements: (i) The amount of charge lost outside a standard extraction 
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aperture, and (ii) the amount of centroid shift experienced by the charge that remains within 
that extraction aperture. This test utilizes the ability of the STIS CCD and its associated 
electronics to read out the image with any amplifier, i.e., by clocking the accumulated charge 
in either direction along both parallel and serial registers. A sequence of nominally identical 
exposures is taken, alternating the readout between amplifiers on opposite sides of the CCD 
(e.g., amps 'B' and 'D' for measuring parallel CTE performance). After correcting for (small) 
gain differences in the two readout amplifier chains, the observed ratio of the fluxes measured 
by the two amps can be fit to a simple CTE model of constant fractional charge loss per 
pixel transfer. Assigning ybin to be the CCD binning factor in the parallel clocking ("y") 
direction and y the pixel location of the target along the y axis where y — 1 at the first pixel 
from the location of amp E0, the measured ratio of signals measured by amps B and D is 



The dependence of the observed signal ratio (e.g., Sb/Sd) on the source position on the CCD 
confirms the consistency with a charge transfer effect (cf. Sections 12.2.11 and 13.1.11 below). 
The CTE is equal to the exponential of half the slope of In (Sb/Sd) vs. y. 

A key virtue of this method is that it requires neither a correction for flat-field response, 
nor an a-priori knowledge of the source flux (as long as the input source is stable during 
the alternating exposures). It should be noted that what is being measured is a sum of the 
charge transfer inefficiencies for the two different clocking directions. However, for identical 
clocking voltages and waveforms and with the expected symmetry of the radiation damage 
effects, we believe the assumption that the CTE loss is equal in the two different clocking 
directions is a reasonable one. 

The implementation of this "internal"^ version of the sparse field test is as follows. Using 
an onboard tungsten lamp, the image of a long and narrow slit is projected at five positions 
along the CCD columns using special commanding of the STIS Mode Select Mechanism. 
At each position, a sequence of exposures is taken, alternating between the 'B' and 'D' 
amplifiers for readout. An illustration of such an exposure sequence is depicted in Figure |2j 
The exposure setup used for these observations is listed in Table [lj while the calibration 
program numbers and dates of each observing epoch are given in Table [2J 



1 This represents the nominal coordinate system of STIS CCD observations in the HST archive 

2 "internal" in this context means that all necessary observations use onboard lamps, so that such ob- 
servations can be performed during Earth occultations, hence not requiring any valuable "external" HST 
observing time 
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The illumination of these images is representative for typical spectroscopic observations 
(as the dispersion direction of STIS CCD spectral modes is essentially along rows). The slit 
image has a narrow profile (2-pixel FWHM), similar to a point source spectrum. The CTE 
loss resulting from this test is "worst-case" , since there is hardly any background intensity 
( "sky" ) to provide filling of charge traps in the CCD array. 

After basic CCD reduction (subtraction of the bias overscan vector and bias image 
subtraction), the average flux per column within a 7-row extraction aperture [which is the 
default extraction size for 1-D spectral extrac t ions i n the CALSTis pipeline used for S TIS 



spectra of point sources, cf. lLeitherer fe Bohlinl Jl997h : lMcGrath. Busko. fe Hodge! Jl999h ] 



as 

well as the centroid of the image profile within those 7 rows are calculated for each exposure. 
In order to reflect the flux measurement method used by the 1-D spectral extraction module 
of CALSTis, a measurement of the background level is obtained 40 pixels above and below the 
extracted flux with a width of 5 rows. Flux and background measurements are performed 
using k — a clipping in order to reject cosmic rays and hot pixels. The alternating exposure 
sequence allows one to separate CTE effects from flux variations produced by warmup of 
the tungsten lamp. As the slit image extends across hundreds of columns, high statistical 
precision on CTE performance can be obtained even at low signal levels per column. 

We emphasize that in calculating CTE loss from this test, charge is only considered 
"lost" if it is no longer within the standard 7-row extraction aperture. I.e., we are only mea- 
suring the component of CTE loss produced by relatively long-time-constant charge trapping. 
Hence, the CTI (= 1 — CTE) values derived from this test will not agree with those measured 
by (e.g.) X-ray stimulation techniques using Fe 55 or Cd 109 , for which charge deferred to even 
the very first trailing pixel contributes to the CTI. However, the measurement described here 
is directly relevant to the estimation of CTE effects on STIS spectrophotometry. 



2.1.2. Spectrophotometric Sensitivity Monitoring 



Another set of data we use to characterize the CTE loss of the STIS CCD in spectro- 
scopic mode is the calibration program that monitors the sensitivity of STIS spectroscopic 
modes. This program consists of wide-slit spectroscopy of primary and secondary spec- 
trophotometric standard stars, using all supported STIS gratings. These data we r e take n on 
a regular basis, with a moni t oring frequency depending on the grating (see iBohlinl (119991 ) and 
Stys. Bohlin. fe Goudfrooijl (120041 ) for details). The particular value of this dataset in the 
context of characterizing the CTE loss is that the spectra of these stars (which are known to 
be intrinsically constant in time) provide a smoothly varying signal level along the dispersion 
direction, thus constraining the signal-level dependence of the CTE loss particularly well. 
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This is discussed further in Section [3T2l 

2.2. Imaging Mode 

2.2.1. External Sparse Field Test 

Similar sparse-field CTE tests using "external" astronomical data have also been carried 
out in Bighf Series of imaging data were acquired on an annual basis for a deld in the 
outskirts of the Galactic globular cluster NGC6752, a field containing several hundreds of 
stars spanning a large range of intrinsic brightness. Every visit of the field consisted of 3 HST 
orbits, in which several exposures were taken using two different exposure times. Several 
repeat exposures were taken at both exposure times, alternating again between opposing 
readout amplifiers. Details of the imaging observations are listed in Table [3J 

We deliberately chose a sparsely populated region in the outskirts of NGC6752, as it 
is well known that the CTE-induced loss in crowded fields is significantly ameliorated (due 
to trap filling) relative to the effects on isolated point sources, while the latter is what we 
attempt to measure here. To allow an assessment of the effect of a varying sky background 
level, we took the data in the so-called Continuous Viewing Zone (CVZ) of HST, in which the 
bright Earth comes closer than usual to the telescope pointing direction. The varying amount 
of scattered light from the bright Earth allows one to obtain a varying "sky" background 
during the CVZ orbits, and hence to obtain measurements of CTI at a suitable range of sky 
background levels. 

Basic CCD reduction steps were performed using tasks in the stsdas.hst_calib.stis 
package of IRAF0. After performing bias and dark subtraction and flatfielding of each indi- 
vidual image using task basic2d, each set of images taken with a given exposure time is 
divided up in subgroups according to the sky background value of the individual images. 
The sky levels are derived by calculating the centroid of the histogram of sky pixel values 
windowed to ±4o" of the mean. This is done using an iterative procedure which rejects 
all pixels outside that window until the number of rejected pixels stays constant from one 
iteration to the next. Each subgroup is subsequently summed together (while performing 
cosmic-ray rejection; task ocrreject), and the summed image is divided by the number 



3 HST Program ID's to date are 8415, 8854, 8911, and 9621. 

4 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Associ- 
ation of Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation, 
U.S.A. 
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of images combined. The sky background values of the resulting "final" images are listed in 
Table EJ 



Aperture photometry is performed using the daophot-ii package (IStetsonl 119871 ) as 
implemented within iraf using fixed-size apertures. Three aperture sizes are used (radii of 
2, 3, and 5 pixels) to evaluate the dependence of calculated CTI values on aperture size (see 
Section [3]). An initial star coordinate list was established using the daofind task, employ- 
ing a detection threshold of 5a above the background. This initial list was then cleaned by 
eliminating targets close to bright stars, in order to avoid anomalous CTI measurements due 
to excessive filling of charge traps. Saturated stars were also eliminated. This cleaned coor- 
dinate list was used for all measurement epochs, using iraf task xyxymatch to calculate 
the appropriate shifts and rotation angles to be applied to the reference coordinate list for 
each epoch. 

Representative results on the parallel CTI for a short-exposure imaging dataset acquired 
in October 2001 are shown in Figure [3] in which the observed flux ratio (amp D/amp B) vs. 
CCD row number is plotted for four different ranges of stellar flux level per exposure. The 
expected CTE behavior is clearly seen, with the closer readout amplifier systematically 
measuring a higher stellar flux than the more distant amplifier. The drawn lines are fits of 
the CTE model we have been considering (i.e., Eq. [Q to the data. Note that the charge loss 
incurred for parallel clocking through the image area of the CCD can be quite substantial. 
For instance, panel (a) in Figure [3] shows that at that time, point sources with a detected 
signal of a few hundred electrons suffer from ~ 25% charge loss when located 1000 rows away 
from the read-out amplifier. Serial CTI values were also determined (by reading out with 
amplifier 'C'), and found to be negligible for all practical purposes (i.e., orders of magnitude 
smaller than the parallel CTI values, and consistent with zero within the uncertainties). The 
default gain = l setting (i.e., 1.0 e~/DN) is used throughout. 

Best-fitting parallel CTI values for different sky background levels are plotted in Fig- 
ure HI This plot illustrates two obvious trends regarding the functional dependence of the 
CTE loss on source signal and background levels: (i) The CTI decreases with increasing 
source signal level, and (ii) At a given source signal level, the CTI decreases with increasing 
sky background levels. These depe ndencies are also seen in investigat ions regarding CTE 



loss of CCDs in the laboratory (e.g.. iHardy. Murowinski. fc Deenlll998l ). 



A comparison of the fits to the data in the different panels of Figure H] shows that the 
slope of log (CTI) vs. log (background) decreases systematically with increasing signal level. 
This suggests (i) that sky background fills traps in the bottoms of the potential wells of the 
CCD, thus mostly ameliorating the transfer of small charge packets, and (ii) a functional 
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dependence of the form 

CTI (B, S) oc exp (-a (B/S) b ) , (2) 

where B is the background level, S is the signal level, and a and b are free parameters. The 
substantial benefit of (only modest) sky background is good news for most science imaging 
observations with the STIS CCD, which typically have longer exposure times than those 
used for these tests and hence will not suffer from the large CTE effects experienced during 
the low-background tests reported here. The functional form of the dependence of CTI on 
signal and background levels will be discussed further in Section [3j 



3. Results 

This section is organized as follows. We first analyse the internal sparse field measure- 
ments to determine the time dependence of CTE loss and to illustrate CTE effects at low 
signal levels when used with the (supported but not commonly used) gain = 4 e~/DN set- 
ting. We then determine the functional forms to correct STIS CCD data for CTE loss as a 
function of observing mode (spectroscopy vs. imaging), signal level, and background level, 
and conclude with an examination of the dependency of the CTE loss on aperture size in 
imaging mode. 



3.1. Internal Sparse Field Data 

3.1.1. The Time Dependence of CTE Loss and Centroid Shift 
The internal sparse field test was the only CTE test conducted with the STIS CCD 



in a uniform manner during both ground testing (at a signal level of ~ 3400 e , iMalumuth 



19961 ) and in-flight operation (on an annual basis). Due to the combination of this long 
baseline in time and the fact that CTI values measured from these data are of great precision 
(as hundreds of columns are averaged together during the analysis), we use these data to 
determine the time dependence of the CTE loss. On the other hand, these data are not 
used to derive the dependences of (spectroscopic) CTE loss on signal and background levels, 
given the fact that these data sample only one (low) background level per signal level, and 
this background level is different for the different signal levels. 

A ^-minimization algorithm is used to compute CTI values for each observing epoch 
and signal level. As outlined in Sect. 12.1.11 the observed ratio of the fluxes measured by 
the two amplifiers was fit to a simple CTE model of constant fractional charge loss per 
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pixel transfer, allowing for k — a clipping of outliers (the latter arise occasionally during the 
short (0.3 s) exposures due to lamp intensity fluctuations). Flux ratio results for the parallel 
internal sparse field test taken after 6.5 years in orbit are presented in Figure [5J As can be 
seen there, the simple CTE model fits the data well. 

To derive the time dependence of the CTE loss, all measurements are first put on a 
uniform scale by normalizing them to one common background level B' = 0. To do so, two 
corrections are required: First, the effect of the spurious charge in STIS CCD bias frames 
(Goudfrooij & Walsh 1997; this is further discussed below) is accounted for by considering 
the total background (£?') to be the measured one (B) plus the spurious charge. Second, 
the background dependence of the CTI according to Eq. [7J (which is discussed in detail 
in Sect. 13.21) is taken into account. The time dependence is then derived by assuming a 
linear increase of CTI with time, as was found earlier for other CCDs aboard HST (e.g., 
Whitmore et al.l Il999l ; lDolphinll2000i ). i.e., by fitting the zero-background CTI values to a 
function of the form: 



CTI(t) = CTI [l + a(* -*„)], (3) 

with t in years and to = 2000.6, the approximate midpoint in time of in-flight STIS observa- 
tions. The conversion between (t — to) and the modified Julian date MJD (which is provided 
by keywords TEXPSTRT and TEXPEND in the science header of each STIS observation) is given 
by t - t = (MJD - 51765)/365.25. 

CTI values derived as mentioned above for the parallel internal sparse field test taken at 
different epochs are plotted in Figure [6j In-flight CTE degradation from a pre-flight starting 
point of low CTI is apparent. Typical CTE behavior is observed as a function of signal level: 
The fractional charge loss (which is proportional to CTI) drops with increasing signal level, 
while the absolute level of charge loss increases. 

Results for the time-dependence fit for gain = 1 are shown in Figure [6] and Table HI The 
functional fit to the data (shown as solid lines) is good, and the derived values for a in Eq. 
2 are consistent with one another (within the uncertainties) for all signal levels measured. 
After weighting the a values of the individual fits by their inverse variance (i.e., RMS" 2 ), the 
weighted average time constant was found to be a = 0.205 ± 0.006, where the uncertainty 
represents the mean error of the mean. 

The effects of charge trapping and release within the 7-row aperture are seen by exam- 
ining the line profiles and centroids. Comparisons of the average line profiles seen for the 
opposing readout directions are shown in Figure [7] for two of the signal levels used. At 60 
electrons per column, the charge trailing and centroid shift are obvious. Even at the higher 



-10- 



signal level, with much lower CTI, the magnified difference between the two observed profiles 
shows that the centroid shift is systematic and measurable. The measured centroid shifts 
(defined as half the difference between the profile centroids as measured by the B vs. the 
D amplifiers) after 5.5 years in orbit are plotted in the left panel of Figure [HI along with 
a least-square fit to the centroid shifts measured at the central position on the CCD as a 
function of measured gross signal level G as read out by the default amplifier D. The right 
panel of Figure [8] depicts the evolution of the centroid shift (at the central location on the 
CCD) due to CTE effects as a function of time. A comparison with Figure shows that the 
centroid shifts increase in time with rates that are quite similar to those measured for the 
CTI values themselves. 

All the above findings indicate that CTE effects cause a fractional loss of signal of 
which the amplitude is directly related to the size of the associated centroid shift, anywhere 
on the CCD. This is illustrated further in Figure [9], which reveals an extremely tight relation 
between measured CTI and centroid shift. The solid line depicts the least-square fit to this 
relation, 

/ CTT \ / CTT \ 2 

Centroid Shift [pixels] = 0.081 — — - 0.002 — — (4) 



Vio- 4 y ' V 10 " 4 / 

which has an RMS error of only 0.01 pixels. (Note that positive centroid shifts in this 
equation indicate shifts toward smaller Y coordinate values in the nominal STIS coordinate 
system.) This relation should prove useful to spectroscopic programs for which both accurate 
flux calibration and astrometry is important. 

Similar measurements have also been carried out in the serial clocking direction by pro- 
jecting vertical slit images at several positions along the CCD rows and alternating between 
amplifiers 'C and 'D' (cf. Figured]) for readout. The resulting charge loss is much smaller 
than for the parallel case and stays below 1% for 1024 pixel transfers, even for the lowest 
signal level measured. This significant difference between the CTE loss experienced in the 
two clocking directions is mainly due to the difference in clocking time, which is a factor ~ 30 
faster in the serial direction f or the STIS CCD (22 /xs pixel -1 in serial vs. 640 /is pixel -1 in 



the parallel direction; see also lBristow &: Alexovll2002l ). Centroid shifts in the serial clocking 



direction are also very small but still well determined, due to the statistical precision offered 
by averaging over hundreds of rows. Interestingly, the serial CTI values and centroid shifts 
are consistent with the relation found above in Eq. 0J This is illustrated in Figure [Ull Given 
the negligible CTE loss experienced by the STIS CCD in the serial clocking direction, the 
remainder of this paper only addresses the parallel CTE loss. 
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3.1.2. CTE Effects at Low Signal Level: Issues for the Gain =4 Setting 

As mentioned above in Section l2.2.1l and shown in Figure HI there is a strong dependence 
of the CTE loss on background level, especially for low signal levels. In fact, this effect is 
significant enough for the structure of STIS CCD bias frames to cause an apparent non- 
linearity in the CTI values derived from the internal sparse field tests at the lowest signal 
levels. This is illustrated in Figure [TT1 The data are shown in two panels per signal level, one 
for each supported CCD gain setting (1.00 and 4.08 e~/DN, hereafter referred to as gain = 1 
and gain = 4 respectively). The open symbols depict actual flux measurements (flux levels 
are indicated on the left ordinate) of the epoch 2002 data of the internal sparse field test. 
The flux measurements with the B and D amps were averaged together for this purpose, after 
having flipped the CCD row numbers for the B-amp measurements around the central row 
(#512) in order to coincide with the D-amp rows. The error bars indicate the measurement 
repeatability. The gain = 4 setting features a CTE-induced loss that is significantly smaller 
than that for the gain = 1 setting. The dashed lines depict the structure of STIS bias frames 
(after subtraction of the overscan vector) in the appropriate gain setting along CCD columns. 
The STIS CCD suffers from accumulative spurious charge, wh ich causes a row-dependen t 



amount of extra "background" that is added during readout (see iGoudfrooij fc Walsh! 1 19971 ). 
The amount of added spurious charge for the gain = 4 setting is much larger than that for 
gain=l (due to a difference in parallel clocking voltages employed), and the slope in the 
spurious charge level across the lower half of the CCD is much steeper for gain = 4 than for 
gain = 1 bias frames. This is the main cause of the difference in low-signal level CTE behavior 
between the two gain settings. Jumping ahead, the filled symbols in Figure [11] depict the 
result of the application of the CTE correction formula (Eq. [7]) derived in Section [3721 below. 
The corrected fluxes are uniform across the CCD to within ~ 1 a for gain = 1, while gain = 4 
data at low signal levels are still somewhat affected by the slope of the spurious charge in the 
bias frames. Hence, the accuracy of the CTE corrections discussed below formally applies 
only to gain = 1 data (and gain = 4 data at signal levels > 2000 e~)@. 



5 The main purpose of the gain = 4 setting of the STIS CCD is to enable very high signal-to-noise obser- 
vations, using signal levels up to the intrinsic full well of the CCD, at the cost of a higher read noise than 
for gain= 1 (see Chapter 7 in lKim Quijano et al.ll2004h . This behavior at low signal levels is therefore not 
detrimental for typical data taken in gain = 4. 
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3.2. Functional Forms for CTE Corrections 

Since spectroscopy and imaging involve inherently different detector illumination pat- 
terns, CTE measurements are expected to show different dependencies on signal level and 
background level for the two observing modes. For instance, a given total signal level within 
a certain aperture size in imaging mode contains a substantial range of signal levels per CCD 
column, and the charge in each of those CCD columns suffers from a different percentage of 
CTE loss. With this in mind, we chose to determine functional forms to characterize the 
CTE loss of the STIS CCD separately for spectroscopy and imaging modes. 



3.2.1. Spectroscopy Mode 

Dependence on Signal and Background Level To constrain the dependence of spec- 
troscopic CTE loss on signal and background levels, we utilize a variety of spectra of spec- 
trophotometric standard stars taken during the last few years of STIS operations, when CTE 
effects were most significant. The standard star spectra used to characterize the CTE effects 
in this report are listed in Table along with their measured signal and background levels. 
All spectra used a 2" wide slit to render slit losses negligible. Spectra are calibrated using 
CALSTis without applying a ny CTE correction, and extracted using traces measured from 



the spectra themselves (see iDressel et al.l ( 120061 ) regarding the evolution of spectral traces 



with STIS). CTI values for spectra of DA white dwarf flux standards GD71 and LDS 749B 
taken using the G230LB grating are derived from the ratio of their measured fluxes by those 
measured from G230L spectra. The (time-dependent) sensitivity calibration for the G230L 
mode (which uses the NUV-MAMA detector of STIS (a microchannel-plate based photon- 
counter, not a CCD) and hence does not suffer from CTE loss) is very well established 



and accurate to subpercent level (jStys et al.l 12004 ). while the sensitivity calibration for the 



CCD spectral modes was established from flux standard star spectra with signal levels high 
enough (> 2 x 10 4 e~ /extraction) to render CTE effects negligible. Figure [121 shows the 
CTI values obtained this way from the G230LB spectra of flux standard stars GD 71 and 
G191— B2B listed in Table |5l which share a common background level. Prior to plotting, 
these CTI values were first normalized to epoch 2000.6 (see Sect. 13.1.11 above). The charge 
loss for low-signal spectra is already quite significant in mid-2000, rising above 10% below 
an extracted signal level of ~ 200 e~ for spectra taken at the nominal central position on the 
CCD. Figure [121 shows that the logarithm of CTI scales roughly linearly with the logarithm 
of signal level, i.e.: 

CTI(G) (xaG~^ (5) 



13 



where G is the (gross) signal level and a and (3 are free parameters. For the G430L and 
G750L spectra in Table CTI values are der ived by comparing the observ e d flux es with pure 



hydrogen white dwarf models, as detailed in lBohlin. Collins, fc Gonnellal (119971 ). The range 



of extracted signal levels covered by the spectra used here is ~ 20 - 16,000 e - pixel -1 , which 
should encompass all spectroscopic observations for which CTE loss is an issue. Smoothed 
versions of the spectra are shown in panel (a) of Figure | 



For the expected basic functional form of the dependence of the CTE loss on signal and 
background levels we combine Eqs. [2] and 

CTI (B, G) = a G- p exp (-7 {B'/G) 5 ) , (6) 

where B' is the total background signal level (i.e., the sum of the sky background B, the dark 
current, and the spurious charge in STIS bias frames which is 0.5 e - pixel -1 for gain = 1 and 
5.0 e" pixel -1 for gain = 4), G is the gross signal level, and a, j3, 7 and 5 are free parameters. 
The values for B and G are re adily obtained from the output of the xlD routine within 
CALSTis to extract 1-D spectra ( McGrath et al.|[l999l ). The background spectrum in xlD is 



smoothed in the spatial and dispersion directions prior to subtraction; this smoothed version 
is the B used in Eq. |6j 

An attempt to fit the function in Eq. O to the CTI values measured for the spectra in 
Table showed that the CTI values of the G230LB and G430L spectra can be fit very well 
this way. However, the CTI values measured for the red end of the wavelength coverage of 
the G750L spectra are systematically found to be below the predictions of this fit. This is 
illustrated in panels (b) and (c) of Figure [T3] which show the best fit of Eq. E]to the G230LB 
and G430L data and the residual percentage of flux error after correcting the data for the 
CTE loss associated with that fit, respectively. In the following we attempt to correct for 
this problem by parametrizing the background parameter in Eq. [6] in terms of the influence 
of the signal in the wings of the point-spread function (PSF) to the filling of charge traps 
during the CCD readout procedure. 



The Impact of the "Red Halo" of the PSF of the STI S CCD The PSF of the STIS 
CCD features broad wings at wavelengths > 8000 A (e.g., iLeitherer fc Bohlinl 119971 ). the 
width of which increases strongly with increasing wavelength. This "red halo" is believed 
to be due to scatter within the CCD mounting substrate which becomes more pronounced 
as the silicon transparency increases at long wavelengths. The effects of the red halo are 
significant, particularly redward of 9 500 A where the default 7 -pixel extraction box captures 
only < 70% of the light in the PSF jLeitherer fc Bohlinl ^997^ . 



This extended halo is likely to have a significant effect on the CTE loss experienced by 
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the signal within the default 7-pixel extraction box, since the charge induced by the halo 
signal that is clocked out before the signal of the source spectrum (i.e., the halo signal on the 
side of the read-out amplifier) should act effectively as 'background' in filling charge traps. 
Note that the red halo signal is not included in the background term B' in Eq. [HI since the 
background spectrum used within calstis/xId is measured far away from the spectrum 
location^]. To improve this situation, we separate the background term into two distinct 
terms, B' (as before) and a new term H which contains the fraction of PSF signal between 
the default 7-pixel extraction box and the read-out amplifier. Values for if as a function of 
wavelength are derived from existing calstis reference files (namely from the Photometric 
Correction Tables^) . We plot these H values in Figure [TH While H is non-negligible at any 
wavelength, the spatial extent of the PSF beyond the default extraction box is only a few 
CCD pixels below ~ 8000 A. Hence, low values of H do not necessarily lower the CTE loss 
significantly. This issue is accounted for in two ways: (%) by subtracting a certain minimum 
threshold value from the measured value of H (i.e., parameter 77 in Eq. [7J below), and (ii) by 
neglecting the halo term for gratings other than G750L or G750M. 

Taking all of the above into account, the full functional form of the CTI algorithm 
considered for spectroscopic modes is 



CTI 5P = aG"" (j(t - 2000.6) + 1) exp ( -5 
where H' = max(0.0, (H — r])) x Net. 



B' + eH' 
G 



(7) 



(Net = G — 7B, the net counts in the extracted spectrum.) Initial estimates of the values 
of parameters a through 77 and their uncertainties were made using bootstrap tests (except 
for the CTI time constant 7, which was fixed at 0.205 as discussed in Sect. 13.1.11 above). 
Particular attention was given to achieving a suitable fit for the full range of signal and 
background levels sampled. A robust fit parameter was then minimized using a non-linear 
minimization routine from Numerical Recipes (Press et al. 1992). Best-fit values of the 
parameters a through 77 are listed in Table [6j 

The quality of this parametrization of the spectroscopic CTE loss is illustrated in Fig- 
ures E] and [T3J The dotted lines in the left panel of Figure E] show the CTI values predicted 
by Eq. [7] for the internal sparse field data at the different gross signal levels (and B' — 0). It 
can be seen that Eq. [7J fits those data within the uncertainties^]. Figure [13] shows the fit of 



6 300 unbimied CCD pixels away by default. 

7 *_pct . f its, listed in data header keyword PCTAB. 

8 The small offset for the 60 e~ /extraction data is due to the effect of the small slope of the spurious charge 
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Eq. [7] to the standard star spectra, both with and without the H term mentioned above. A 
comparison of panels (c) and (e) of Figure [13] in particular clarifies that the inclusion of the 
H term yields a significantly better correction for the red end of G750L spectra. Overall, 
the new parametrization of the CTE loss is accurate within 5% for any data point in the 
spectra used here, while the RMS accuracy for all spectra used in this study together stays 
within 1%. For reference, the dotted lines in panel (e) in Figure [TBI depict the Poisson noise 
associated with a resolution element (assumed to be 2 pixels along the dispersion) in spectra 
of a given signal level. This CTE correction formula renders spectroscopic flux calibration to 
an accuracy better than the uncertainty due to Poisson noise per resolution element. 

3.2.2. Imaging Mode 

CTI values are calculated for every epoch of the external sparse field test as explained in 
Sect. 12.2. 1| grouped by ranges of flux and sky background levels. The ranges of signal levels 
are chosen after inspecting histograms of signal levels for each sky level listed in Table [3J 
aiming for similar numbers of targets in each group. Table [7] lists the derived CTI values for 
each group and epoch. 

Time Dependence: Consistency With Spectroscopic Mode Before embarking on 
a comprehensive procedure to fit the imaging CTI values as a function of signal and back- 
ground levels, we derive the time dependence of the CTE loss from these data and check for 
consistency with the time constant derived before in Sect. 13.1.11 for spectroscopic mode. To 
this end, we select groups of signal and sky levels that were measured consistently in at least 
three epochs and whose CTI uncertainties were smaller than the CTI values themselves. The 
CTI values of these groups are plotted in Figure [TS] as a function of time, along with least- 
squares fits of CTI vs. time (separately to each group). After weighting the slopes of the 
individual fits by their inverse variances, the weighted average time constant is 0.214 ± 0.047, 
i.e., consistent with the 0.205 ± 0.006 found in Sect. 13.1.11 We adopt the latter in the follow- 
ing, given the much longer time baseline and the more homogeneous and well-determined 
data used to derive it. 

Dependence on Signal and Background Level The functional form to characterize 
the CTI of the STIS CCD in imaging mode was determined in a way similar to that of the 



level in the lower half of STIS CCD bias frames, as discussed in Sect. 13.1.21 and illustrated in Figure [TT] 
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spectroscopic mode. We first attempted to fit the imaging CTI values with the function 
in Eq. [6] used to fit the spectroscopic CTI values, but found that the dependence on back- 
ground level required more terms in this case. After extensive experimentation, the following 
functional form produced the best fit to the imaging data: 

CTI/m = a e~ b lcts x (c (t - t ) + 1) x 

where 'counts' stands for the net counts within the measurement aperture, and 

t-t = (MJD — 51765)/365.25 

lcts = ln(counts) — 8.5 

bck = max(0.0, sky) 

lbck = ln(sqrt(bck 2 + 1)) - 2 

The constants 51765, 8.5, and 2 were roughly the averages of the corresponding parameters 
in the data, and were included to provide numerical stability as well as to produce indepen- 
dent coefficients (a through g). The purpose of the 'bck' and 'lbck' parameters is to avoid 
logarithms of negative values. (The parameter c was fixed to 0.205 as discussed above.) The 
fitting procedure was the same as that used in Sect. 13. 27T1 above. The resulting best-fit values 
of the coefficients in Eq. [8] are listed in Table [8J 

The quality of this parametrization of the CTE correction is depicted in Figure [16j 
Quantitatively, it yields a correction that is accurate within 4<r of any single data point in 
Tabled while the overall RMS uncertainty is 3.0%, which is equivalent to 0.2 a of the "mean" 
data points plotted in panel (b) of Figure dHJ To put this in perspective, an observation of a 
typical faint star in September 2002 with an extracted signal of 100 e~ and a background of 
6 e~ per pixel at the center of the CCD underwent a CTE loss of ~ 15% (cf. Figure [T5|) . The 
parametrization in Eq. [8] corrects this loss to an RMS photometric accuracy of only 0.5%, 
much lower than that due to Poisson noise for that object. 



( ( bck y\ 

dexp (— e lbck) + (1 — d) exp I — / I 

\ V counts/ / 



Centroid Shifts in Imaging Mode Centroid shifts were also measured for the external 
sparse field imaging data, defined (as before) as half the difference between the centroids of 
the target stars as measured by the two different amplifiers. 

The derived centroid shifts are plotted against measured CTI in Figure [T7J which shows 
a strong relation between the two parameters, similar to the case for the spectroscopic data 
described in Sect. 13.1.11 The solid line depicts a weighted least-square fit to this relation for 
the imaging data, 

Centroid Shift [pixels] = 0.025 ( J - 0.78 10" 3 I — 
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which has an RMS error of 0.005 pixelg^]. We emphasize that centroid shifts due to CTE 
effects are rarely taken into account (a litera ture search using the NASA Astrophysics 
Data System^] came up with only one study: Piatek et al. 2005 ). whereas they are rele- 
vant for science programs that aim to reach the highest possible positional accuracy such as 
proper motion measurements for which accuracy of order 0.01 pixel is often required (e.g., 



McNamara. Harrison, fc Andersonll2003l ; iKoslowski et al.ll2006l ). This is especially important 



if one combines measurements done with two different instruments (or with one instrument 
but at significantly different locations on the CCD) to derive proper motions. 

A comparison of Eq. [9] with Eq. [U(or, equivalently, Figure [T71 with Figure Mj shows that 
the imaging data incur a roughly 3 times smaller centroid shift than spectroscopic data for a 
given CTI value. We interpret this difference in the sense that the measured centroid shifts 
are most heavily influenced by the centroid shift incurred by the peak pixel of the PSF (since 
it is by far the brightest pixel within the measurement aperture), and higher signal levels 
cause lower centroid shifts. To give an illustrative estimate of this effect: A comparison of 
Eq. [8] with Eq. [7] shows that for a given CTI value and background level, the signal level in 
imaging mode is a factor ~ 7 higher than in spectroscopic mode. Approximating the imaging 
and spectroscopic PSFs by Gaussians with a FWHM of 1.8 pixels, this translates to a ratio 
in peak intensity 

PeakiMA 7 

PeaksPEC ^/^/{A\n2) FWHM PSF 

According to Eq. a 3.6 times higher signal level causes a ~ 3.0 times lower CTI value, and 
hence a 2.9 times lower centroid shift (cf. Eq. HJ), which is consistent with what we find in 
practice. 



The Effect of Aperture Size for Imaging Mode Since the amplitude of the centroid 
shift increases with CTE loss as discussed above, one obvious question is whether this effect 
has an impact on small-aperture photometry, which is often used in imaging programs of 
compact sources in order to achieve an increased signal-to-noise ratio relative to the use of 
larger apertures. We checked for this by evaluating the dependence of the measured CTE 
loss on aperture size using the imaging data of the external sparse field test. Photometry 
measurements were made through three aperture sizes (radii of 2, 3, and 5 pixels). The 
result is illustrated in Figure [18j There is no significant difference in CTI value among the 
three apertures used: The formal ratio between CTI values measured through 3-pixel radius 



'Again, positive centroid shifts in this equation indicate shifts toward smaller Y coordinate values, 
'http: / /www. adsabs.harvard.edu 
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and through 2-pixel radius apertures for the data shown in Figure [TBI is 0.998 ± 0.015 (mean 
error of the mean). Likewise, the ratio between CTI values measured through 5-pixel radius 
and through 2-pixel radius apertures is 1.031 ± 0.020. We interpret this result in the sense 
that the centroid shift to the PSF caused by CTE effects does not cause actual loss of signal 
as long as one measures the centroids of the targets on the image itself with an adequately 
large centering box widtli"]. 



3.3. Implementation of CTE Correction 

As we have adopted the definition of CTE as the fraction of charge transferred per pixel 
transfer, the correction for CTI (= 1 — CTE) to be applied to an extracted STIS data point 
OBJ (be it a flux element of a STIS CCD spectrum 12 ! or an entry of a star in a photometry 
table from STIS CCD imaging) is 

OBJcorr = — j- (10) 

(1 - CTIobj) 1 ^-y*yhn 

where y is the Y coordinate of the centroid of the target, ybin is the on-chip binning factor, 
and CTIobj is the result of Eq. [7| [for spectroscopic data] or Eq. [H] [for imaging data] for that 
data point. 



4. Concluding Remarks 

Using a variety of methods, a set of formulae to correct observations of point sources 
with the STIS CCD for CTE-related loss of signal has been derived and presented. For 
data taken in imaging mode, the CTE loss is parametrized in terms of the location of the 
source on the CCD, the source signal level within the measurement aperture, the background 
level, and the time of observation. For spectroscopic data, one extra parameter is needed 
to provide an adequate calibration of the CTE loss, namely the signal in the point spread 
function located between the signal extraction box and the read-out amplifier. The effect of 
the latter parameter is significant at wavelengths redward of 8000 A where the PSF of the 
STIS CCD exhibits a extended "halo", i.e., for spectra taken using the G750L or G750M 
gratings of STIS. The algorithms presented here correct flux calibration inaccuracies due to 



11 We used the default size of the centering box within daophot-ii (cbox = 5 pixels). 

12 Note however that since December 16, 2003, the pipeline flux calibration of STIS CC D spectra takes 
correction for CTI already into account for full- frame datasets ( Goudfrooii fc Bohlin 20061 ) . 
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CTE losses as large as 30% to within ~ 1.5% RMS throughout the wavelength range covered 
by the STIS CCD observing modes. 

Using bi-directional CCD readouts, centroid shifts incurred due to CTE loss are also 
derived. Centroid shifts are up to ~ 1CT 1 pixel for imaging photometry of faint point sources, 
and a factor of ~ 3 higher for spectroscopy mode. A strong correlation is found between the 
signal loss and the centroid shift (both for imaging and spectroscopic modes), thus enabling 
one to correct for both effects of imperfect CTE to STIS CCD observations. While centroid 
shifts due to CTE effects are rarely taken into account in the literature, they are relevant for 
science programs that aim to reach the highest possible positional accuracy such as proper 
motion measurements for which accuracy of order 0.01 pixel is often required. 

Though we have demonstrated that the correction formulae work well, it is of course 
preferable (when practical) to minimize the magnitude of the CTE effects on an observing 
program at the data acquisition stage. To this end, STIS operations support acquisition 
of targets for CCD spectroscopic modes at the so-called 'El pseudo-apertures', which are 
located ~ 125 rows from the readout amplifier rather than at field center. The appropriate 
instrumental calibrations have been developed for these aperture positions. At the STIS 
plate scale of 0.05 arcseconds per pixel along the slit, this position gives 6 arcseconds of sky 
coverage to the near edge of the CCD; targets for which this angular coverage is sufficient 
can thus be observed with four times lower charge loss by selecting this observing position 
rather than the center of the CCD. Corrections for CTE loss and the use of the 'El pseudo- 
apertures' will be of even greater importance in the future if STIS is successfully repaired as 
planned in the hoped-for HST Servicing Mission 4, as the CCD usage will be extended for a 
number of additional years in the orbital radiation environment. 
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Fig. 1. — Schematic architecture of the STIS CCD. The 1024 x 1024 pixel device has two 
serial registers and four readout amplifiers. The image can be read out through any single 
amplifier or via two- and four-amplifier combinations. The nominal amplifier (amp D) is at 
the top right. 



Fig. 2. — Representative images used for the "internal" sparse field CTE test in the parallel 
clocking direction. At each of the five positions along the CCD columns, a sequence of 
exposures is taken, alternating between amplifiers 'B' and 'D' (see Figured]). Systematic 
variation of the relative signals measured by the two amplifiers as a function of position 
reveals the CTE effects. The images are displayed at a logarithmic stretch. 



-24- 



— i 1 1 1 1 1 1 1 1 r— i 1 1 1 1 1 1 1 1 

14 b- 200 < flux < 500, background 3 e 

0.6 ™ T (a) 
—1 — 1 — 1 — I — 1 — 1 — 1 — I — 1 — 1 — 1 — I — 1 — 1 — 1 — I — 1 — 1 — 1 — \ 



2 1 4 - 600 < flux < 1500, background 3 e 
R * h -4 



< 1.2 



1 

0.8 
0.6 



CTI = (1.4 ± 0.1) 10 

(b) 

— 1 — 1 — 1 — I — 1 — 1 — 1 — I — 1 — 1 — 1 — I — 1 — 1 — 1 — I — 1 — 1 — 1 — \ 



1.2 

PQ 
(X 

a 0.8 

\ 0.6 

O 

g 1,4 £- 6 .° 00 < flux < 15,000, background 3 e 
CTI = (5.2 ± 0.3) 10" 5 



(c)1 

_i 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 L 



200 400 600 800 1000 

Row nr. from B amp 



Fig. 3. — Flux ratio results for the external imaging sparse field CTE test in the parallel 
clocking direction, taken in October 2001 (after 4.7 years in orbit). The sky background 
was 3 e~ /pixel. The three panels correspond to measurements at different signal levels, as 
shown in the legends. In each panel, the ratio of the signals measured by the two different 
amplifiers is plotted as a function of the Y position of the image along the CCD columns. 
The best-fitting simple CTE model (constant fractional charge loss per pixel transfer) as 
well as the corresponding CTI value are shown in each panel. The stellar fluxes shown here 
were measured using a circular aperture with a radius of 3 pixels. 
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Fig. 4. — Log-log plot of CTI vs. sky background level for six different signal levels, derived 
from the external imaging sparse field test taken after 4.7 years in orbit. The flux levels, the 
slopes of the best-fitting line, and the latter's uncertainties (in parentheses) are mentioned 
in each panel. Note the systematic decrease of the CTI dependence on the sky background 
with increasing signal level. 
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Fig. 5. — Flux ratio results for the 2002 epoch of the internal sparse field CTE test in the 
parallel clocking direction. The six panels correspond to measurements at different signal 
levels, as indicated above and below each panel. In each panel, the ratio of the signals 
measured by the two different amplifiers is plotted as a function of the mean Y position of 
the image along the CCD columns. The best-fitting simple CTE model (constant fractional 
charge loss per pixel transfer) as well as the corresponding CTI value are shown in each panel. 
Star symbols indicate measurements used in the fit and circles indicate rejected points. 
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Fig. 6. — Left panel: CTI normalized to zero background for gain = 1 as a function of time 
and signal level, derived from the internal sparse field test. Symbols associated with the data 
of individual signal levels (corrected for CTI) are indicated in the legend. The solid lines 
represent linear fits to the data (one for each signal level), while the dotted lines represent the 
result of applying Eq. [7] to the data (one line per signal level group). Right panel: Absolute 
charge lost due to CTI for an object at the central row of the STIS CCD as a function of 
time and signal level. Symbol types are the same as in the left panel. The epoch of HST 
Servicing Mission 2 (during which STIS was installed on HST) is depicted as a black dotted 
vertical line. 
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Fig. 7. — Comparison of line profiles obtained when clocking a given exposure in opposite 
directions in the internal sparse field test. At low signal levels (left panel), charge trailing 
and centroid shift are obvious. Even at high signal levels (right panel), differences in the 
line profile are systematic and measureable (the dotted line depicts the difference profile 
multiplied by fO). Figure reproduced from Kimble, Goudfrooij, & Gilliland (2000). 
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Fig. 8. — Left panel: The centroid shift (in unbinned CCD pixels) as a function of signal level 
as read out by the D amplifier for the gain = 1 observing block of the internal sparse field test 
in October 2002, ~5.5 years after STIS installation. Centroid shifts for the central location 
on the CCD are shown in filled squares, and a least-squares fit to the latter is shown by the 
solid line. Right panel: The centroid shift for the central location on the CCD as functions 
of time and signal level as read out by the D amplifier. Both the data and the corresponding 
linear fits are plotted. Symbols associated with individual signal levels (corrected for CTI) 
are indicated in the legend. 
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Fig. 9. — The centroid shift at the central row of the CCD as a function of CTI for all 
gain = l datasets of the internal sparse field test. Symbols as in Figure [HI Note the very 
tight relation between the two properties. A least-squares fit to the data (cf. Eq. HJ) is shown 
by the dashed line. 
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Fig. 10. — Panel (a): Serial CTI as a function of signal level for the 1999 epoch of the 
internal sparse field test. Panel (b): The centroid shift at the central column of the CCD as 
a function of CTI for the same data. The relation between CTI and centroid shift found for 
the parallel CTI (Eq. H]), as shown in Figure [9]) is shown by the dashed line. 
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Fig. 11. — Results of the two lowest-signal cases of the 2002 epoch data of the internal sparse 
field test in the parallel clocking direction, for the two gain settings (1 and 4 e~/DN). The 
average signals measured at the five different positions along the CCD columns are shown 
as symbols (open circles for gain = 1 (left panels), open squares for gain = 4 (right panels)). 
Extracted signal levels are shown on the left-hand ordinate. The dashed lines depict the 
'background' level due to spurious charge in STIS CCD superbias frames, averaged over all 
columns. Superbias intensity levels (in e~ per pixel) are shown on the right-hand ordinates. 
Note the strong effect of local 'background' to the CTI for these low signal levels, especially 
for gain = 4 data. The latter show a much lower CTI than the gain = 1 data, due to a 
significantly higher amount of spurious charge in gain = 4. The filled symbols depict the 
signal levels after correction for CTI using Eqs. [7] and [TUJ 
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Electrons in 7-pixel extraction 

Fig. 12. — Fractional charge loss at the central row of the CCD (i.e., (1 — CTI) 512 ) vs. 
extracted gross signal level for three G230LB wide-slit spectra of flux standard stars observed 
at the center of the STIS CCD. Crosses, squares, and triangles represent values for spectra 
o8v204030, 08v2040e0, and o6igl0010, respectively (see Table [5]). The CTI values were 
normalized to correspond to epoch 2000.6 (see Sect. 13.1 .ID . Note the linear dependence of 
log (CTI) on log (signal level). 
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Fig. 13. — Panel (a): Smoothed flux standard star spectra used to determine the functional 
form of the CTI of the STIS CCD in spectroscopic mode. The legend links the symbols with 
the STIS rootnames of the individual spectra, as listed in Table |5j Panel (b): Fractional 
charge loss at the central row of the CCD vs. gross signal level within the default 7-pixel 
extraction box. Symbols represent measured CTI values for the spectra shown in panel (a), 
normalized to epoch 2000.6 (see Sect. 13.1.11) . The smooth curves represent the predictions 
of the best-fit CTI model without 'red halo' term [i.e., Eq. [6] for those data. Symbol types 
are the same as in panel (a). Panel (c): The ratio of measured CTE values and the model 
predictions shown in panel (b) vs. gross signal level. Panel (d): Same as panel (b), but 
now using the final CTI model including a 'red halo' term [i.e., Eq. [7J. Panel (e): Same 
as panel (c), but now using the final CTI model including a 'red halo' term. For reference, 
the dashed lines represent the uncertainty due to Poisson noise associated with the binned 
spectra shown in panel (a) as a function of signal level, while the dotted lines represent the 
Poisson error associated with a resolution element (2 pixels along the dispersion) of unbinned 
spectra. 
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Fig. 14. — Parameter H in Eq. [7J The fraction of the light in the PSF above the default 
7-pixel extraction box as a function of wavelength for the STIS CCD grating modes. Note 
the discontinuity near 3000 A, at the boundary of the wavelength ranges covered by the 
G230LB and G430L gratings. This is likely due to the presence of a Lyot stop in the G430L 
and G 750L modes which is absent in the G230LB mode (see also lDressel et al.ll2006l ; [Proffitt 
2006h . 
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Fig. 15. — Imaging CTI as a function of time, derived from the external sparse field test. 
Both the data and the corresponding linear fits are plotted (the latter in dotted lines). 
Symbols associated with individual signal levels are indicated in the legend (at the top left). 
Data with different sky background have symbols with different colors: Black for sky = 3 e~ 
pixel -1 , red for sky = 6 e" pixel -1 , blue for sky = 10 e - pixel -1 , and magenta for sky = 30 
e - pixel -1 . 
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Fig. 16. — CTI at epoch 2000.6 as a function of parameter 'lets' (as defined in Sect. 13.2.21) 
for imaging data derived from the external sparse field test. Panel (a) shows all individual 
data points listed in Table El while panel (b) shows weighted mean data points (averaged 
over 5 ranges of 'lets', namely —3 < lets, —3 < lets < —2, —2 < lets < —1, —1 < lets < 0.5, 
and lets > 0.5, separately for each sky level). The drawn lines in both panels depict the 
fit of Eq. [8] to the data. Symbol and line colors associated with individual sky levels (in e~ 
pixel -1 ) are indicated in the legend in panel (a) (at the bottom left). 
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Fig. 17. — The centroid shift at the central row of the CCD as a function of CTI for the 
imaging data taken during the external sparse field tests. A weighted least-squares fit to the 
data (cf. Eq. [9]) is shown by the solid line. 
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Fig. 18. — The ratio between imaging CTI values measured in different aperture sizes. 
Left panel: 3-pixel radius vs. 2-pixel radius. Right panel: 5-pixel radius vs. 2-pixel radius. 
Open circles indicate measurements from the Sep 2002 epoch, while filled squares indicate 
measurements from the Sep 2003 epoch. The size of the measurement aperture does not 
have any significant impact on the CTI in imaging mode. 
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Table 1. STIS observations for "internal" sparse field tests 



Optical 


■^ccn 


Slit 


Exp. Time 


CCD Gain 


Signal 


Element 


(A) 


Name 


(s) 


(e-/DN) 


Levels a 


G430M 


5471 


0.05x31NDA 


0.3 


1,4 


60 


G430M 


5471 


0.05x31NDA 


0.6 


1,4 


130 


G430M 


5471 


0.05x31NDA 


0.9 


1,4 


195 


G430M 


5471 


0.05x31NDA 


2.3 


1,4 


500 


Mirror 


b 


0.05x31NDA 


0.3 


1,4 


3400 


Mirror 


b 


0.05x31NDB 


0.3 


1,4 


9800 



a Unit is e per column in the absence of CTE losses. 
b Using "Clear" filter (i.e., no filter; aperture name 50CCD) 

Note. — All exposures were taken in two read modes: Once read 
out using the default amplifier "D", and once using amplifier "B", 
on the opposite side of the parallel register. 

Table 2. Observing blocks used for the internal sparse field test. Each block extended 
over a time period of one to a few days. Representative values for the Modified Julian Date 

(MJD) and the civil date are shown. 



Gain 


Block 


Program 


MJD 


Date (UT) 


1 


1 


8414 


51428 


Sep 07, 1999 


1 


2 


8414 


51650 


Apr 16, 2000 


1 


3 


8851 


51845 


Oct 28, 2000 


1 


4 


8910 


52210 


Oct 28, 2001 


1 


5 


9620 


52567 


Oct 20, 2002 


1 


6 


9620 


52896 


Sep 14, 2003 


4 


1 


8414 


51783 


Aug 26, 2000 


4 


2 


8851 


52004 


Apr 04, 2001 


4 


3 


8910 


52399 


May 05, 2002 


4 


4 


9620 


52734 


Apr 05, 2003 


4 


5 


10026 


53099 


Apr 04, 2004 
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Table 3. STIS observations 51 for "external" sparse imaging field test 



Observing 


On-orbit Time 


Exp. Time 


N cxp 


CCD Gain 


Background 


Date 


(yr) 


(s) 




(e~/DN) 


Levels b 


Sep 15, 1999 


2.54 


20 


24 


1 


3.1, 5.1 






100 


9 


1 


14 


Oct 13, 2000 


3.62 


20 


40 


1 


3.3, 5.9, 7.9 






100 


12 


1 


16, 27 


Sep 13, 2001 


4.54 


20 


40 


1 


3.0, 5.8, 7.9 






100 


12 


1 


11, 37 


Sep 13, 2002 


5.45 


20 


40 


1 


2.9, 5.3 






100 


12 


1 


11, 33 


Sep 13, 2003 


6.51 


20 


40 


1 


3.4 






100 


12 


1 


11, 16 



a The imaging field in NGC 6752 was centered on a = 19 d ll m 05 ! ;96, 5 = 
— 60°01' 06'.'7 in J2000 equinox. Observations were conducted using the "Clear" 
filter (i.e., no filter; aperture name 50CCD). 

b Centroids of histogram of sky pixel values windowed to ±4cr of the mean. 
Unit is e~ per pixel. 



Note. - - All exposures were taken in two read modes: Once read out using 
the default amplifier "D" , and once using amplifier "B" , on the opposite side of 
the parallel register. Celestial coordinates are given in equinox J2000. 

Table 4: CTI time constant a as measured from the internal sparse field test data taken in 
gain = 1, as a function of signal level. 



signal 


a 




(DN) 


(y^ 1 ) 


(yr- 1 ) 


60 


0.216 


0.009 


130 


0.192 


0.013 


195 


0.188 


0.021 


500 


0.202 


0.016 


3450 


0.218 


0.038 


9850 


0.170 


0.052 
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Table 5: List of flux standard star spectra used to characterize the CTI as function of signal 
and background level. All intensities are in e~ per pixel per (sub-) exposure. 



Rootname 


Grating 


Flux Standard 


Observing 


Background 


Range in 








Date 


level 


Signal Levels 


o6igl0010 


G230LB 


G191-B2B 


2002-04-22 


0.4 


1000 - 5000 


o6ig01060 


G750L 


GD71 


2001-12-07 


1.7 


300 - 16000 


o6il01020 


G230LB 


LDS 749B 


2001-07-13 


1.9 


100 - 1800 


o8u2200t>0 


G430L 


AGK+81D266 


2003-11-03 


0.5 


3000 - 9200 


o8vl01030 


G750L 


WD 1657+343 


2004-04-07 


2.5 


30 - 750 


o8v2040e0 


G230LB 


GD71 


2003-10-24 


0.3 


120 - 730 


o8v204030 


G230LB 


GD71 


2003-10-24 


0.1 


20 - 170 



-43- 



Table 6: Best-fit Values of Coefficients in Functional Form of Spectroscopic CTI (Eq. [7]). 



Coefficient 


Value 


Description 


a 


0.056 ± 0.001 


CTI normalization 


P 


0.82 ± 0.01 


Gross count level dependence 


7 


0.205 ± 0.006 


Time dependence of CTI 


<5 


3.00 ± 0.05 


Multiplicative factor for 'background'/gross count ratio 


e 


1.30 ± 0.10 


Multiplicative factor for halo light fraction 


c 


0.18 ± 0.01 


Power of 'background'/gross count ratio 


7] 


0.06 ± 0.01 


Minimum value of halo light fraction above spectrum 
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Table 7: Measured CTI values and centroid shifts for external sparse imaging field test data. 



MJD 


Sky Level 
(e~ / pixel) 


Signal Level 
(e") 


CTI 

(pixel - 


1 x 

x ) 




Centroid Shift 
(pixel) 


51436 


3.1 


149 ± 27 


(2.3 ± 


0.4) 


10" 


-4 


0.0793 


± 


0.0268 




3.1 


293 ± 87 


(1.8 ± 


0.2) 


10" 


-4 


0.0485 


± 


0.0124 




3.1 


460 ± 82 


(1.9 ± 


0.2) 


10" 


-4 


0.0397 


± 


0.0114 




3.1 


852 ± 276 


(1.4 ± 


0.1) 


10" 


-4 


0.0359 


± 


0.0084 




3.1 


2293 ± 567 


(8.3 ± 


0.5) 


10" 


-5 


0.0226 


± 


0.0056 




3.1 


2888 ± 593 


(7.6 ± 


0.5) 


10" 


-5 


0.0211 


± 


0.0063 




3.1 


4708 ± 913 


(5.7 ± 


0.4) 


10" 


-5 


0.0183 


± 


0.0070 




5.1 


290 ± 85 


(1.4 ± 


0.2) 


10" 


-4 


0.0418 


± 


0.0110 




5.1 


476 ± 84 


(1.5 ± 


0.2) 


10" 


-4 


0.0318 


± 


0.0115 




5.1 


836 ± 274 


(1.2 ± 


0.1) 


10" 


-4 


0.0268 


± 


0.0073 




5.1 


2320 ± 569 


(7.2 ± 


0.4) 


10" 


-5 


0.0192 


± 


0.0055 




5.1 


2902 ± 587 


(6.5 ± 


0.4) 


10" 


-5 


0.0182 


± 


0.0053 




5.1 


4765 ± 876 


(5.4 ± 


0.4) 


10" 


-5 


0.0158 


± 


0.0060 




14.1 


915 ±279 


(5.7 ± 


0.8) 


10" 


-5 


0.0225 


± 


0.0068 




14.1 


2231 ±532 


(6.3 ± 


0.6) 


10" 


-5 


0.0100 


± 


0.0048 




14.1 


4833 ± 838 


(3.9 ± 


0.4) 


10" 


-5 


0.0111 


± 


0.0042 




14.1 


10721 ±2157 


(3.4 ± 


0.3) 


10" 


-5 


0.0086 


± 


0.0028 




14.1 


23397 ±8367 


(2.5 ± 


0.2) 


10" 


-5 


0.0055 


± 


0.0019 


51831 


3.3 


413 ±101 


(2.1 ± 


0.2) 


10" 


-4 


0.0436 


± 


0.0157 




3.3 


951 ± 193 


(1.4 ± 


0.1) 


10" 


-4 


0.0276 


± 


0.0069 




3.3 


1141 ±308 


(1.5 ± 


0.1) 


10" 


-4 


0.0267 


± 


0.0073 




3.3 


3055 ± 659 


(9.9 ± 


0.5) 


10" 


-5 


0.0181 


± 


0.0033 




3.3 


3586 ± 689 


(9.1 ± 


0.4) 


10" 


-5 


0.0161 


± 


0.0035 




3.3 


5600 ± 957 


(7.9 ± 


0.4) 


10" 


-5 


0.0120 


± 


0.0034 




3.3 


11594 ±2812 


(5.2 ± 


0.2) 


10" 


-5 


0.0086 


± 


0.0023 




3.3 


26075 ± 8100 


(3.4 ± 


0.2) 


10" 


-5 


0.0072 


± 


0.0028 




5.9 


407 ±119 


(1.4 ± 


0.2) 


10" 


-4 


0.0340 


± 


0.0158 




5.9 


931 ± 195 


(1.1 ± 


0.1) 


10" 


-4 


0.0309 


± 


0.0084 




5.9 


1156 ±315 


(1.1 ± 


0.1) 


10" 


-4 


0.0265 


± 


0.0076 




5.9 


3049 ± 675 


(8.0 ± 


0.4) 


10" 


-5 


0.0197 


± 


0.0050 




5.9 


3661 ± 689 


(7.1 ± 


0.4) 


10" 


-5 


0.0179 


± 


0.0053 




5.9 


5585 ± 925 


(6.8 ± 


0.3) 


10" 


-5 


0.0169 


± 


0.0052 




5.9 


11666 ±2821 


(4.8 ± 


0.2) 


10" 


-5 


0.0116 


± 


0.0048 
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Table 7: (continued) 



MJD Sky Level Signal Level 



CTI 



Centroid Shift 



(e /pixel) 


(e ) 


(pixel 


) 




(pixel) 


7.9 


417±113 


(1.4 ± 


0.2) 


10" 


-4 


0.0315 




0.0128 


7.9 


930 ± 183 


(1.0 ± 


0.1) 


10" 


-4 


0228 




0.0076 


7.9 


1153 ±303 


(1.1 ± 


0.1) 


10" 


-4 


0.0224 




0.0064 


7.9 


3007 ±661 


(7.6 ± 


0.4) 


io- 


-5 


0.0199 




0.0034 


7.9 


3622 ±691 


(7.1 ± 


0.4) 


10" 


-5 


0.0180 




0032 


7.9 


5546 ± 945 


(6 5 + 


31 


io- 


-5 


n 01 38 


-|- 




7.9 


11656 ±2819 


(4.4 ± 


0.2) 


10" 


-5 


0.0104 


± 


0.0019 


7.9 


25093 ± 6925 


(3.3 ± 


0.1) 


10" 


-5 


0.0067 


± 


0.0020 


14.8 


1188 ±325 


(9.2 ± 


0.1) 


io- 


-5 


0.0199 


± 


0.0067 


14.8 


2846 ± 718 


(6.8 ± 


0.8) 


10~ 


-5 


0.0130 


± 


0.0042 


14.8 


3604 ± 742 


(5.8 ± 


0.9) 


io- 


-5 


0.0144 


± 


0.0046 


14.8 


6047 ± 937 


(5.3 ± 


0.4) 


io- 


-5 


0.0139 


± 


0.0037 


14.8 


13100 ±2780 


(3.9 ± 


0.3) 


io- 


-5 


0.0103 


± 


0.0029 


14.8 


30780 ±11882 


(2.9 ± 


0.1) 


10" 


-5 


0.0052 


± 


0.0017 


29.5 


1197 ±344 


(6.1 ± 


0.8) 


10" 


-5 


0.0169 


± 


0.0058 


29.5 


2809 ± 714 


(5.9 ± 


0.6) 


10" 


-5 


0.0137 


± 


0.0040 


29.5 


3403 ± 696 


(4.9 ± 


0.8) 


10" 


-5 


0.0124 


± 


0.0049 


29.5 


6268 ±917 


(5.1 ± 


0.4) 


10" 


-5 


0.0091 


± 


0.0038 


29.5 


13111 ±2866 


(3.3 ± 


0.3) 


10" 


-5 


0.0087 




0.0018 


29.5 


30806 ±11776 


(2.5 ± 


0.1) 


10" 


-5 


0.0054 


± 


0.0012 


3.0 


148 ± 30 


(3.0 ± 


0.2) 


10" 


-4 


0.0831 


± 


0.0178 


3.0 


302 ± 88 


(2.4 ± 


0.2) 


io- 


-4 


0.0604 


± 


0.0102 


3.0 


904 ± 253 


(1.4 ± 


0.1) 


10~ 


-4 


0.0371 


± 


0.0069 


3.0 


3031 ±882 


(1.0 ± 


0.1) 


10" 


-4 


0.0251 


± 


0.0036 


3.0 


4739 ± 1598 


(7.7 ± 


0.5) 


10" 


-5 


0.0224 


± 


0.0036 


3.0 


10969 ± 3939 


(5.2 ± 


0.3) 


10" 


-5 


0.0147 


± 


0.0049 


3.0 


29354 ±11702 


(4.4 ± 


0.3) 


10" 


-5 


0.0091 


± 


0.0036 


5.8 


149 ± 30 


(2.1 ± 


0.3) 


io- 


-4 


0.0486 


± 


0.0197 


5.8 


281 ±84 


(1.8 ± 


0.2) 


10" 


-4 


0.0387 


± 


0.0108 


5.8 


893 ± 254 


(1.2 ± 


0.1) 


10" 


-4 


0.0238 


± 


0.0066 


5.8 


3046 ± 901 


(8.8 ± 


0.4) 


io- 


-5 


0.0202 


± 


0.0043 


5.8 


4811 ±1612 


(7.1 ± 


0.5) 


10" 


-5 


0.0197 


± 


0.0060 



51831 



52166 
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Table 7: (continued) 



MJD 


Sky Level 


Signal Level 




CTI 






Centroid Shift 




(e — /pixel) 


(e") 


(pixel 


) 




(pixel) 


52166 


7.9 


149 ±31 


(2.1 




0.2) 


io- 


-4 


0.0378 




0.0182 




7.9 


292 ± 83 


(1.5 


± 


0.1) 


10" 


-4 


0.0341 




0.0100 




7.9 


912 ±255 


(1.1 


± 


0.1) 


10" 


-4 


0.0248 




0.0054 




7.9 


3042 ± 901 


(8.6 


± 


0.3) 


10" 


-5 


0.0197 


± 


0029 




7.9 


4812 ±1620 


(7.1 


± 


0.4) 


10" 


-5 


0.0188 


± 


0033 




7.9 


10801 ± 4001 


(4.9 


± 


0.3) 


10" 


-5 


0.0103 


± 


0.0039 




7.9 


29533 ±11383 


(4.3 


± 


0.3) 


10" 


-5 


0.0050 


± 


0.0023 




11.4 


385 ± 76 


(1.3 


± 


0.2) 


10" 


-4 


0.0287 


± 


0.0141 




11.4 


977 ±236 


(1.0 


± 


0.1) 


10" 


-4 


0.0230 


± 


0.0060 




11.4 


3217 ±867 


(5.5 


± 


0.5) 


10" 


-o 


0.0130 


± 


0.0032 




11.4 


4818 ±1126 


(4.9 


± 


0.3) 


10" 


-o 


0.0098 


± 


0.0032 




11.4 


10145 ±3215 


(4.3 


± 


0.2) 


10" 


-u 


0.0085 


± 


0.0020 




11.4 


35503 ± 23516 


(2.6 


± 


0.2) 


10" 


r: 

-o 


0.0054 


± 


0.0019 




35.2 


960 ± 235 


(6.1 


± 


0.8) 


10" 


c 

-o 


0.0093 


± 


0.0044 




35.2 


3225 ± 905 


(4.5 


± 


0.4) 


10" 


K 

o 


0.0084 


± 


0.0036 




35.2 


4864 ± 1139 


(4.2 


± 


0.4) 


10" 


K 

o 


0.0077 


± 


0.0026 




35.2 


10218 ±3289 


(3.9 


± 


0.2) 


10" 


r: 
-u 


0.0064 


± 


0.0025 




35.2 


34799 ± 23530 


(2.7 


± 


0.1) 


10" 


c 

-o 


0.0042 


± 


0.0017 


52499 


2.9 


127 ±35 


(A 7 


-|- 


41 


10" 


-4 


fl D823 


-|- 


fl D247 




2.9 


182 ±37 


f3 5 


-|- 


n 3l 


10" 


-4 


fl D721 


-|- 


n m 53 




2.9 


369 ± 103 


(3.4 


± 


0.2) 


10" 


-4 


0.0584 


± 


0.0105 




2.9 


1131 ±321 


(1.8 


± 


0.1) 


10" 


-4 


0.0276 


± 


0.0045 




2.9 


3833 ± 1076 


(1.2 


± 


0.5) 


10" 


-4 


0.0198 


± 


0.0034 




2.9 


5610 ±1198 


(9.4 


± 


0.6) 


10" 


-5 


0.0178 


± 


0.0036 




2.9 


12893 ±4821 


(6.0 


± 


0.4) 


10" 


-5 


0.0140 


± 


0.0056 




5.3 


185 ±37 


(2.8 


± 


0.4) 


10" 


-4 


0.0913 


± 


0.0213 




5.3 


362 ± 108 


(2.6 


± 


0.2) 


10" 


-4 


0.0504 


± 


0.0122 




5.3 


461 ± 76 


(2.0 


± 


0.2) 


10" 


-4 


0.0488 


± 


0.0124 




5.3 


1082 ± 297 


(1.5 


± 


0.1) 


10" 


-4 


0.0394 


± 


0.0055 




11.0 


520 ±132 


(1.6 


± 


0.2) 


10" 


-4 


0.0349 


± 


0.0101 




11.0 


1226 ± 305 


(1.1 


± 


0.1) 


10" 


-4 


0.0298 


± 


0.0054 




11.0 


4029 ± 1095 


(6.8 


± 


0.5) 


io- 


-5 


0.0225 


± 


0.0031 




11.0 


6037 ±1438 


(6.2 


± 


0.4) 


10" 


-5 


0.0165 


± 


0.0029 




11.0 


13224 ±4429 


(5.4 


± 


0.3) 


10" 


-5 


0.0128 


± 


0.0025 
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Table 7: (continued) 



MJD Sky Level Signal Level 


CTI 


Centroid Shift 


(e - /pixel) (e - ) 


(pixel -1 ) 


(pixel) 



52499 



52885 



32.6 


483 ±111 


(1 


1 


± 





3) 


10" 


-4 





0310 


± 





0126 


32.6 


507 ±133 


(9 


1 


± 





2) 


10" 


-5 





0332 


± 





0108 


32.6 


1237 ±312 


(6 


2 


± 





8) 


10' 


-5 





0251 


± 





0057 


32.6 


4067 ±1108 


(5 


2 


± 





5) 


10" 


-5 





0239 


± 





0033 


32.6 


5986 ± 1344 


(5 


1 


± 





4) 


10" 


-5 





0173 


± 





0035 


32.6 


13158 ±4413 


(4 


6 


± 





2) 


10" 


-5 





0145 


± 





0022 


32.6 


42713 ±8317 


(2 


6 


± 





2) 


10" 


-5 





0104 


± 





0014 


3.4 


141 ± 27 


(4 


o 


± 


o 


4) 


10" 


-4 


o 


0979 


± 


o 


0308 


3.4 


250 ± 55 


(3 


5 


± 


o 


3) 


10" 


-4 


o 


0863 


± 


o 


0152 


3.4 


372 ± 107 


(2 


9 


± 


o 


2) 


10" 


-4 


o 


0709 


± 


o 


0125 


3.4 


468 ± 75 


(2 


3 


± 


o 


3) 


10" 


-4 


o 


0643 


± 


o 


0153 


3.4 


1118 ±310 


(1 


9 


± 


o 


1) 


10" 


-4 


o 


0473 


± 


o 


0058 


3.4 


3862 ± 1061 


(1 


3 


± 





1) 


10" 


-4 





0329 


± 





0043 


3.4 


5630 ± 1207 


(1 





± 





1) 


10' 


-4 





0270 


± 





0050 


3.4 


14176 ±4217 


(5 


6 


± 





5) 


10" 


-5 





0145 


± 





0038 


3.4 


34584 ± 13298 


(4 





± 





3) 


10" 


-5 





0085 


± 





0031 


11.3 


517 ±74 


(2 





± 





2) 


10" 


-4 





0442 


± 





0130 


11.3 


1212 ±305 


(1 


3 


± 





1) 


10" 


-4 





0294 


± 





0046 


11.3 


4077 ±1100 


(8 


1 


± 





5) 


10" 


-5 





0236 


± 





0033 


11.3 


6155 ± 1432 


(7 


4 


± 





4) 


10" 


-5 





0204 


± 





0034 


11.3 


15001 ±4020 


(5 


6 


± 





2) 


10' 


-5 





0134 


± 





0019 


11.3 


35832 ± 18767 


(3 


3 


± 





2) 


10' 


-5 





0076 


± 





0015 


15.5 


540 ± 137 


(1 


4 


± 





3) 


10' 


-4 





0337 


± 





0153 


15.5 


1213 ±304 


(1 





± 





1) 


10" 


-4 





0295 


± 





0045 


15.5 


4064 ± 1084 


(6 


3 


± 





5) 


10" 


-5 





0211 


± 





0039 


15.5 


6086 ± 1416 


(6 


3 


± 





4) 


10" 


-5 





0176 


± 





0036 


15.5 


15225 ±4272 


(5 


4 


± 





3) 


10" 


-5 





0139 


± 





0021 


15.5 


37163 ±18658 


(3 


1 


± 





2) 


10" 


-5 





0082 


± 





0022 
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Table 8: Best-fit Values of Coefficients in Functional Form of Imaging CTI (Eq. [H]). 



Coefficient 


Value 


Description 


a 


(1.33 ± 0.03) 10~ 4 


CTI normalization 


b 


0.54 ± 0.01 


Count level dependence 


c 


0.205 ± 0.006 


Time dependence of CTI 


d 


0.05 ± 0.02 


Fraction of dependence on background (vs. background/flux ratio) 


e 


0.82 ± 0.05 


Multiplicative factor for background level dependence 


f 


3.60 ± 0.07 


Multiplicative factor of background/flux count ratio dependence 


9 


0.21 ± 0.02 


Power of background /flux count ratio 



